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Abstract. We conjecture an integrability and linearizability test for dispersive Z^-lattice equa- 
tions by using a discrete multiscale analysis. The lowest order secularity conditions from the 
multiscale expansion give a partial differential equation of the form of the nonlinear Schrodinger 
(NLS) equation. If the starting lattice equation is integrable then the resulting NLS equation 
turns out to be integrable, while if the starting equation is linearizable we get a linear Schrodinger 
equation. On the other hand, if we start with a non-integrable lattice equation we may obtain 
a non-integrable NLS equation. This conjecture is confirmed by many examples. 

1. Introduction 

By a difference equation we mean a functional relation, linear or nonlinear, between functions 
calculated at different points of a lattice [4,5,11,23,39]. These systems appear in many appli- 
cations. First of all they can be written down as discretizations of a differential equation when 
one is trying to solve it with a computer. In such a case one reduces the differential equation to 
a recurrence relation: 

du 

— = f{x,u), xGM =^ v{n + 1) = g{n,v{n)), n G Z. 
ax 

On the other hand we can consider dynamical systems defined on a lattice, i.e. systems where 
the real independent fields depend on a set of independent variables which vary partly on the 
integers and partly on the reals. For example we can consider the differential-difference equation: 

(fu{n, t) 



F{t, u{n, t),u{n — l,t), ..,u{n — a, t),u{n + 1, t), .., u{n + b, t)), 



with n, a, 6 S N, t G M. These kind of equations can appear in many different settings. Among 
them they are associated with the evolution of many body problems, the study of crystals, 
biological and economical systems and so on. 
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The world of the differential-difference equations is indeed much richer than that of partial 
differential equations as one can see in [32], where it is shown that almost any Hamiltonian 
network of weakly coupled oscillators has a "breather" solution, while the existence of breathers 
for a nonlinear wave equation is rare. These results imply that the discrete world can be richer 
of interesting solutions and thus worthwhile studying by itself. 

In the discrete world one can easily use computers to solve the relevant discrete maps but 
their solution is always approximate due to rounding numerical errors. So the possibility of 
treating discrete models by exact techniques is very important as in this way new phenomena 
may be uncovered. Exact techniques are naturally associated with integrable models, so that to 
find techniques for showing the integrability of discrete systems is a challenging task. 

Various different criteria have been proposed as tests for integrability in the discrete setting. 
One of the earliest proposals for ordinary difference equations was the singularity confinement 
test of Grammaticos, Ramani and Papageorgiou [12], which has proved to be an extremely 
useful tool for isolating discrete Painleve equations. Hietarinta and Viallet, after discovering 
some non-integrable equations with the singularity confinement property, were led to introduce 
the algebraic entropy condition for integrability of rational maps [17] (the phenomenon of weak 
degree growth of integrable maps had been studied earlier by Veselov [40,41] and it is connected 
with the Arnold complexity [1]). Ablowitz, Halburd and Herbst extended the Painleve property 
to difference equations using Nevanlinna theory [3] by considering the asymptotic growth of 
meromorphic solutions at infinity. Roberts and Vivaldi [35] have studied the distribution of 
orbit lengths in rational maps reduced to finite fields ¥p for different primes p, in order to 
identify integrable cases of such maps. The algebraic entropy has been the only method applied 
to detect integrability also in the case of two-dimensional partial difference equations (or Z^- 
lattice equations) [18,43]. 

In the continuous setting multiscale techniques [37, 38] have proved to be important tools 
for finding approximate solutions to many physical problems by reducing a given dispersive 
nonlinear partial differential equation to a simpler equation, which is often integrable [7]. These 
multiscale expansions are structurally strong and can be applied to both integrable and non- 
integrable systems. Zakharov and Kuznetsov in the introduction of their article [45] say: " If 
the initial system is not integrable, the result can be both integrable and nonintegrable. But 
if we treat the integrable system properly, we again must get from it an integrable system". 
Calogero and Eckhaus [7] used similar ideas starting from a class of hyperbolic systems to prove 
in 1987 the necessary conditions for the integrability of dispersive nonlinear partial differential 
equations. Later Degasperis and Procesi [10] introduced the notion of asymptotic integrability 
of order n by requiring that the multiscale expansion be verified up to a fixed order n. 

Recently a few attempts to carry over this approach to difference-difference and differential- 
difference equations have been proposed [6,13,14,20,24-26,28]. In [20,25,28] we developed 
a multiscale expansion technique on the lattice which, starting from dispersive integrable Z^- 
lattice equations, provided non-integrable Z^-lattice equations, thus contradicting the Zakharov 
and Kuznetsov's claim. 

Later on, in [13], this problem has been solved by extending the previous results to func- 
tions of infinite slow-varyness order (see Section [2] for details). This new technique easily fits 
for both difference-difference and differential-difference equations. The representative example 
considered in [13] was the lattice potential KdV equation. In such a case a proper representa- 
tion of the discrete shift operators in terms of differential operators, which is equivalent to a 
Taylor expansion, provides an integrable nonlinear Schrodinger (NLS) equation as the lowest 
order secularity conditions from the multiscale expansion. Results obtained by performing the 
Taylor expansion of the shifted variables can be found in many papers, see for instance [8,22]. 
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However, it is worthwhile to note that our method always allows us to get (for free) a discrete 
equation from a continuous one by going over to a finite slow-varyness order. 

The result found in [13] suggests the investigation of a discrete analogue of the Calogero- 
Eckhaus theorem [7] , which claims that a necessary condition for the integrability of a dispersive 
nonlinear partial differential equation, is that its lowest order of the multiscale expansion be an 
integrable differential equation. 

In the present work, by using the formalism developed in [13], we propose the following 
conjecture: 

Conjecture. If a nonlinear dispersive 7? -lattice equation is integrable (resp. lineariz- 
ahle) then its lowest order multiscale reduction will he an integrable nonlinear (resp. linear) 
Schrddinger equation. 

To confirm the above Conjecture we shall consider several examples, both integrable and non- 
integrable, linearizable and non-linearizable. In Section[2]we present a brief review of the discrete 
multiscale method, developed in [13]. Then, in Section [3] we apply the multiscale expansion to 
some difference-difference KdV-type and Toda-type equations. In Section S] we consider some 
linearizable equations, as the differential-difference Burgers equation and the Hietarinta equa- 
tion, and also a non-linearizable difference-difference Burgers-type equation. Finally, Section [5] 
is devoted to open problems and concluding remarks. 

2. Basic formulas for the discrete multiscale analysis 

The basic tool of the discrete reductive perturbation technique developed in [13] is a proper 
multiscale expansion which consists in considering various lattices and functions defined on them. 
The relation between the variation of a function on two different lattices of indices n and ni is 
given by [21] 

Kun = ^(-ir^ P n„+, =i! j;^Ai^n„„ (1) 

1=0 ^ ^ i=j 

where u„ : Z — > M is a function defined on a lattice of index n ^ 7L and : Z — s- M is 
the same function on a lattice of index n\ G Z. According to Eq. ([T]), by the symbol A„ 
we mean the standard forward difference of the function Un with respect to its subscript, e.g. 

Un — lin+i ~ Un, where Tn is the shift operator TnUn — Un+i- The coefficients 
Pij in Eq. ([T]) are expressed in terms of the ratio of the lattice spacing for the variable ni with 
respect to that of variable n. 

Eq. (pP) implies that a finite difference in the discrete variable n depends on an infinite number 
of differences on the variable ni, e.g. the function Un+i can be written as a combination of the 
functions itj's, for i varying on an infinite set of points of the lattice ni. In [13] one has considered 
generalizations of formula ([T]) in order to deal with functions n„ = u^.^^.^k ^ depending on a 
finite number K of lattice variables Ui and with functions depending on two discrete indices, say 
n and m, thus dealing with Z^-lattice equations. 

To get a reduction of a given difference equation onto a difference equation of order less than 
a fixed number, say i, one has to consider functions Un of slow-varyness of order i, namely the 
space of those functions Un such that A^^Un = or, equivalently, Af^^Um = (see [28]). With 
this definition one can reduce the infinite series expansion ([T]) to a finite number of terms. 

To deal with functions of infinite order of slow-varyness one considers a formal expansion of 
the shift operator T„. By introducing on the lattice of index n the real variable x = nax, the 
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shift operator such that T^u^x) = u{x + a) can be formally written as 



Introducing a formal derivative with respect to the index n, say 6n, we can define the discrete 
operator Tn as 

oo 

Tn = eMSn)=Y,^. (2) 
The formal expansion ([2]) can be inverted, yielding 

Sn = lnT„, = ln(l + A„,) = V ^^A^, (3) 

i=l 

where A„ is the discrete first right difference operator with respect to the variable n, see Eq. 
([T]). We refer to [13] for more technical details and decompositions with respect to different 
discrete derivatives. 

To perform a multiscale expansion we need to consider functions defined on different lattices, 
thus depending on a fast lattice index n and on multiple slow- varying lattice indices nj, 1 < 
i < K. The slow- varying lattice variables vary on a larger scale with respect to the one of the 
original lattice of index n, and thus the transition from n to ni corresponds to a coarse graining 
of the lattice. In the continuous limit, when the spacing between the lattice points goes to zero, 
this corresponds to the introduction of multiple continuous variables: given x € M we define the 
new variables Xi = e*x, 0<e<^l, l<i<K. By taking into account the above definitions 
we can introduce a function u r \k = u{x; {xjj^i ) depending on a fast index n and K slow 

indices rii = e*n. Here e = 1/N and N is an integer number if we require that the slow indices 
have to be integer numbers. 

At the continuous level, the total derivative dx acting on functions u{x;{xi}^i) is the sum 
of partial derivatives, i.e. dx = dx + edx^ + e^S^j + 0{e^). Consequently we can expand the total 
shift operator Tx in terms of the partial shift operators 

Tx = ex.p{axdx) = exp{axdx) ex.p{eaxdx) eyip{e'^axdx2) ■■■ ■ 



At the discrete level, we can write 

Tn = exp((5„) exp(e(5„i) exp(e^(5„2) • • • = TnT^^Tf^ 



with 

oo oo ^ oo 22 

T — \ " IIL — \ La* T^^ — \" ^ 

i=0 ■ i=0 ' 1=0 

where the (5-operators are given in Eq. ([3]) and they act on a function u^.^^.^k ^ with respect 
to their subscript. From the previous formulas we deduce immediately that working with 5- 
operators is equivalent, up to terms depending on ax, to perform a Taylor expansion of the 
discrete function ^i„.{„.}if ^, e.g. for = 2 we have 



g2 

1 + e6n, + -Sl^ + e^6n2 + O(e^) 



In the following Sections we shall discuss a list of examples of integrable, non-integrable, 
linearizable and non-linearizable dispersive lattice equations. 
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3. MULTISCALE EXPANSION OF DISCRETE TODA-TYPE AND KdV-TYPE EQUATIONS 

3.1. An integrable difference-difference Toda equation. An integrable discrete-time Toda 
equation is given by [27]: 

exp(Un,m - Un,m+l) " exp(n„,m+i - Un,m+2) = (4) 
= a^[exp(-Un-l,m+2 - Un,m+l) " exp(-Un,m+l " Un+l,m)], 

where a G R is a parameter related to discretization of the time variable. The above five-point 
Z^-lattice equation has been obtained by Hirota thirty years ago [19]. 

We can split Eq. (jlj) into a linear and nonlinear part by considering its small amplitude 
solutions, namely Un^m = ^Wn,m, < e <C 1, e being a small parameter as the one we considered 
in Section [TJ The linear part of Eq. ^ has a travelling wave solution of the form Wn,m = 
exp {i[Kn — a;(«;)m]}. Here the dispersion relation lo = lo{k) obeys to the following equation: 

a2(f7-/C)2 = /C(17-l)2, (5) 

where we have introduced the quantities K, = exp(iK) and Vt = exp(— ia;), G M. 

As nonlinearity generates harmonics, we introduce the following expansion for the function 



W„ 



aeZ k=0 



where wj^ = w^. , w being the complex conjugate of w. Moreover, to avoid secularities we 

have to require that w^^^ = w^j^\ni, {mi}fLi), where ni and {mi}fLi are K + 1 slow-varying 
lattice variables, namely ni = en and rrii = e^m, 1 < i < K . 

The multiscale expansion of Eq. (jH) provides several determining equations for the coefficients 
obtained selecting the different powers of e and the different harmonics a. Let us give here 
the main results obtained by considering the lowest e and a orders of this development. 

At 0(1), for a = 0, 1, we find linear equations which are identically satisfied either directly 
or by taking into account the dispersion relation ([5]). For |a| > 2 one gets some linear equations 

whose unique solution is given by tt^g"^ = 0. 

At 0{e), for the harmonics a = 1,2, we find the following equations: 

Mni + <5^J4'^ = 0, (7) 

-P=rirfV, (8) 

where 

{n-i){n + }C) ic + i 

'^0 = ^rTTTP^ — — ; '^1 ^ 



217(/C-1) ' ' 2(/C-l)' 
By solving Eq. @ with respect to a; = ui{k) one gets that tq = Vg = duj/dn, where Vg is the 
group velocity. The solution to Eq. ([7]) is given by w^Q\ni,{mi}f^i) = w^Q\n2,{'mi}f^2) with 
722 = ni — Vgirii. Eq. ([8]) expresses w^^^ in terms of w^^\ Moreover, as is a function of n2 

the same must hold for w^f^ , i.e. w'^^\ni, {mi)f^^) = w'f\n2, {"ii}^2)- 
We can now consider the 0{e^). For a = we have: 

[(r2C +^3(C -25n/n^J]u;^ =T4(5„, -25„J|4^)|2, (9) 

with 

6^^(/C-l)(l^^-/C) Qjcnjn - If 2 

= ' = " [^-icy ' ^4 = 6(0-1). 
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As w^^ is a function of n2, the same must be for w'^q\ Thus we can integrate Eq. ([9]) by 
requiring that its solution be bounded and we get 

[{T2V^g + T3(l + 2Vg)]6n,wS'^ = ^{1 + 2Vg)\w^o^\^ (10) 

For a = 1 we find a secular equation for which is solved by requiring that satisfies 
the following equation: 

fx I x2 ^ (1) _ (0) I -(1) (2) n-w 

\T5Om2 + T6f„2jW^5 = T7W^ 0n2WQ + TsWq ' W{ ' , [LI) 

where 

12(/C- 1)^2 3(172 _/C)(0- 1)2 



T5 = 7^ z: , T6 



T-j = 4tq, Ts 



n-jc ' " 2(/c-i) 
6{ic + - jc){n - 1)2 



Taking into account Eqs. (|8I1U|) . Eq. (fTT|) reduces to the following NLS equation for w^^: 



where 

sin(K + u) + sin a; — sin(2u; + k) 



T"9 



no 



8(cos K — \) 

sin(2u; + k)(cos k + 5) — 10sin(w + /t/2) cos(k/2) 



4(cos K — 1) 

As Tg and tiq are real parameters the NLS equation (|12p is integrable and it can be obtained 
as a compatibility conditions for a Lax pair, see [14]. 

As noticed in [13], assuming a finite order of slow-varyness in the NLS equation (|12|) we get 
a non-integrable difference-difference equation. In other words, a necessary condition for the 
integrability of Eq. (jl2|) is that d. = oo. 

3.2. A non-integrable difference-difference Toda equation. The multiscale expansion of 
the standard Toda lattice, 

iln = exp(n„_i - Un) - exp(u„ - n„+i), (13) 

has been carried out by Kalyakin in [22], after the publication of the pioneering work by Zakharov 
and Kuznetsov [45]. Here Un = Un{t), and by dot we mean the time derivative. 

We shall perform the multiscale analysis of a non-integrable dispersive discretization of 
Eq. (US). It reads 

Un,m+1 ^Uj-i^m ~l~ '^n/m— 1 — 0,[eXJp(^'Uj-i~l,m ^7i,m) ^^p{,Ufi,m 'Wn+l,m)]i (1^) 

where a G M is a parameter related to the discretization of the time variable. 

We shall proceed, as in the previous example, by considering Eq. (jl4p for small amplitudes 
Un,m = GU!n,m- In this casc the linear part of the resulting equation admits the dispersion relation 

a;(K) = 2 arccos asin^— ^ . (15) 

Using the expansion ([6]) we get a set of determing equations similar to the one we obtained 
for the integrable difference-difference Toda equation (jlj), but with different coefficients. 
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The resulting NLS equation for the harmonic 'wl^\ obtained at 0{e^), is again given by 



Eq. p2p with the following real coefficients: 



Vg cos u) — a cos K 
rg = - 



no 



2 sin w 
a(cos K — 1) 
sinw 



2a(cos K — 1) sin^ k 

+ cos K - 1 + 



v"^ — a (a — f)(cos k — 1) 



where Vg = dto/dK is the group velocity corresponding to the dispersion relation (|15p . 

3.3. Non-integrable difference-difference KdV equations. We now present the multiscale 
analysis of two different non-integrable lattice KdV equations obtained by discretizing the con- 
tinuous KdV equation. They read: 

'Un,m+l 'Un,m— 1 — ^ ('Un+3,m ^'^n+l,m ~l" ^Un—\^rn "Wn— 3,m) (-^^6) 



h 

--[u, 



and 



2 

a 



2 

n+l,m 



/2 ^ 



_^ / 2 _ 2 \ 

where a, b are real parameters. Clearly, the difference between Eqs. (|16p and (|17p lies just in 
their nonlinear part. Precisely, Eq. (|16p has a symmetric nonlinear part, while it is asymmetric 
in Eq. pzp . Since the linear part of both Eqs. (|16p and (jl7p is the same, they have the same 
dispersion relation, given by 

uj{k) = arcsin(asin^ k). 

Let us start with the analysis of Eq. (jl6p by looking at small amplitudes solutions Un^m = 
ewn,m and considering the expansion ([6|). 

At 0(1) we get a set of determining equations which give w^"^ = for \a\ > 1, while the 
0{e) gives w^"* = and the following deter ming equations: 

Mm =0' (18) 
ii;f)=ai(4^))2, (19) 

where 

du) b cos K 

^5 = ^1 = — . 2 q \ • 

clK iasm k(4cos'^ K — coswj 

As in the previous examples, Eq. (fTH]) is solved by w'o\ni, {mi}f^^) = w^\n2, {mi}^2) with 
712 = rii — VgUii. Moreover, higher harmonics imply that = if |a| > 3. 



At C(e^) we get 



4°)=a2rf^P, a2 = -, (20) 



for a = 0. For a = 1, after removing secularities and using Eqs. ()19|20p we find the following 
NLS equation: 

iSm2Wo^ = (^s^l^w^o^ + (^iw^o^ \wo^ l"^ , (21) 
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where 



<73 



asmK[3(l — 3cos^ k) — — cos^ k)] 



(74 



a cos Lo sm k 



2 cos a; 

cos OJ COS n 

+ 



SCOSK 2(cos u — 4 COS'^ K 



It is remarkable to notice that Eq. (j2ip is an integrable NLS equation (if d. = oo) even if the 
starting difference-difference KdV equation (jl6p is non-integrable. At this level the integrability 
of Eq. (j2ip is definitely due to the symmetric version of the nonlinear part of Eq. ()16p . As a 
matter of fact, the same multiscale expansion carried out on Eq. (jl7p yields a NLS equation of 
the form (j2ip . but with complex coefficients 0"3,cj4, thus breaking its integrability property. In 
this latter case the coefficients Uj, 1 < i < 4, read 



and 



6exp(iK;) h 

^1 = ~A ^ Ta ^ ^2 = T, ! 

4asin k(4cos'^ K — cosa;j 2,Vg 
a sin k[3 - 9 cos^ - vi s\v? k] \h{ai + (72) [1 - exp(i«;)] 

cos LO I COS LO 

4. Multiscale expansion of linearizable discrete equations 

In this Section we shall consider difference-difference and differential-difference equations, 
which are linearized by a transformation of coordinates, usually a contact transformation. In 
particular, we analyze a differential-difference Burgers equation and the Hietarinta equation. Let 
us here recall that a multiscale reduction of the Hietarinta equation has been already considered 
in [28] for £ = 2, while its linearizability has been proven in [34]. 

4.1. Discrete Burgers equations. We look for a discrete Burgers equation which has a real 
dispersion relation. This request is equivalent to say that its linear part admits dispersive wave 
solutions. A discrete Burgers equation has been derived in [15], but it is not dispersive since 
it involves an asymmetric discrete-time derivative. One can prove that a dispersive completely 
discrete Burgers- type equation which is linearizable via a discrete Cole-Hopf map does not exist. 
Roughly speaking this is related to the fact that integrability implies non-symmetric discrete- 
time derivatives, which are usually incompatible with real dispersion relations. 

Therefore we shall shall consider the following dispersive differential-difference Burgers equa- 
tion [30]: 

ia iln = [1 + aUn){Un+l - Un) + — , (22) 

1 + aUn-1 

with a G M. Eq. ()22p is obtained as a compatibility condition for the following Lax pair: 

a{l + aun-i) 

As Eq. (I22p has a real dispersion relation we can use the same procedure described in Section 
[H but with a continuous-time variable. The linear part of Eq. (I22p admits wave solutions with 
a dispersion relation given by 

U!(k) = — ^(COS K — 1). 

For solutions of Eq. (f^^ with small complex amplitude M„(t) = ewn{t), where Wn{t) is 

expanded according to a natural modification of Eq. ([6]), we get at 0(1) the conditions w^^^ = 
for a > 1 and a < 0. 
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At 0{e) we have = for a > 2 and a < and w^^'^ = 0. Moreover we get 

{vg6m + 5ti)w^^^ = 0, 

whose solutfon is wl^\ni, {ti}^^) = wl^\n2, {tj}^2) with n2 = ni — Vgti, Vg = duj/dn, and 

= (24) 



2 cos K — COS 



At 0(6^) we get w^"^ = for a = 0, while for a = 1, using Eq. (j24|) . we obtain a hnear 
equation 

8„„,a) = (25) 

a remainder of the fact that Eq. (j22p is hnearizable. 

Let us now consider the multiscale analysis of a difference-difference equation obtained 
through a naive, but symmetric, discretization of Eq. (I22p . Obviously the resulting equation is 
not hnearizable and it is not associated with a Lax pair of the type given in Eq. (I23p . It reads 

• 2 

7^(Wn,m+l - ■"n,m-l) = (1 + au„ (^in+l,m " + —T-f — ■, (26) 

lb 1 + aUn-l,m 

where 6 € M is a parameter related to the time discretization. The linear part of Eq. (j26p admits 
dispersive wave solutions with dispersion relation 



uj{k) = — arcsin 



2b 

^(COS K 



By applying the usual procedure, we get, at O(e^), the analogue of Eq. ([251) . It reads 

(1) _ a^Vg[Vg{cos K - 1) -cosk] (i) 

which is again a linear equation. 

4.2. The Hietarinta equation. The Z^-lattice equation 

,m+l + 02 

has been introduced by Hietarinta in [16]. Here, Cj and Oj, i = 1,2, are real and distinct 
parameters. As proven in [34], Eq. ()27p is hnearizable. For our purposes it is convenient to 
reparametrize the coefficients Cj and Oj by means of the following transformations: Cj i— > and 
Oi 1-^ o~^, 1 = 1,2. 

As noticed in [28] the linear part of Eq. (I27p has travelling wave solutions Un,m = exp{i[Kn — 
Lo{K,)mW with a complex dispersion relation. Nevertheless, the constraint oi + ei = 02 + 62 
provides the real dispersion relation 



(27) 



u:{k) = 2arctan 



a — b ( K 
tan 



a + b \2 



where a = 02 — oi, 6 = ei — 62. 

For small amplitude solutions Un,m = £'Wn,m, with Wn,m given by Eq. ([6]) we obtain from 
Eq. ([27|) a set of equations which, at 0(1), give the conditions zug"^ = for a > 1 and a < 0. 

At C(e) we find wf'' = and 

[VgSn, + dm,)w\'^ = 0, (28) 
«^?^=Ci(4'V, (29) 
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where 



du> exp(iK) + exp[i(K + uj) 

Vg = -r-. Ci = ei + o- 



dn exp(iK) — exp(ia;) 

Eq. (f28]l is solved by w^i\ni, {ti}f^i) = w^^\n2, {ti}f=2) with n2 = rei — Vgti. 
At order ©(e^), by taking into account Eq. (f29]) . we get for a = 0, 

wf = (2ei + 6-a)|4^)|2, 
which leads, for a = 1, to the following linear equation 

i6m2Wi^ = Cs ^n2'^i\ Cs = ^(cos K - COS sincj. (30) 

In our previous work [28], the multiscale analysis of the Hietarinta equation provided, in the 
i = 2 case, the following difference-difference equation: 

i('/'n2,m2+l 4'n2,m2^ ~l~ Cl ((/'rt2+2,m2 20^2, m,2 ~l~ 4'n2—2,m2 ) + 

n2+l,m2 2(f)n2,m2 + 4'n2 — l,m2 ) + C3(j)n2 ,m2 \ 4'n2,m2 r + (31) 
n2,m2) 4^n2,m2 — 0; 

where (j) = w^'^ and the function ip is defined through the equation 

V'n2+l,m2 - n2 — l,m2 — [(/)ji2,m2 (0n2+l,™'2 4'n2 — l,m2) ~^ 4'n2,m2i4'n2+l,m2 4'n2 — l,m2 

)]. (32) 

The coefficients Cj, 1 < i < 6, have been determined and they are all real. Eq. ()3ip is a (non- 
integrable) discrete integral equation due to the presence of the function ■0. It is remarkable 
to note that, if £ = oo, Eq. ([32]) can be explicitly integrated and the resulting coefficients q 
combine in such a way that the nonlinear term in Eq. (I3ip vanishes, thus giving Eq. (I30p . 



5. Concluding remarks 

In this paper we presented a Conjecture which gives necessary conditions for the integrability 
and linearizability of dispersive difference-difference and differential-difference equations. This 
Conjecture has been confirmed by a long list of examples, contained in Sections [3l [H In par- 
ticular, performing the discrete multiscale analysis, we found that the resulting NLS equation 
turns out to be integrable if: 

(1) the starting Z^-lattice equation is integrable, or 

(2) the starting Z^-lattice equation is non-integrable, but obtained through a symmetric 
discretization of a continuous integrable equation. 

We expect that claim (2) fails at higher orders of the multiscale expansion. Moreover the re- 
sulting NLS equation turns out to be non-integrable if the starting nonlinear Z^-lattice equation 
is obtained through a non-symmetric discretization of a continuous integrable equation. Let us 
recall that discrete symmetric derivatives play a fundamental role in the integrability of discrete 
equations, see [44]. We stress here that the integrability of the NLS equation is lost as soon as 
we go to a finite slow-varyness order since shifts operators do not satisfy the Leibniz rule and 
thus the reduced Lax pair is not compatible anymore. 

Concerning the linearizability of Z^-lattice equations we obtained a linear Schrodinger equa- 
tion if: 

(1) the starting Z^-lattice equation is linearizable, or 

(2) the starting Z^-lattice equation is obtained through a discretization of a continuous 
linearizable equation. 
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The above characterization is evident just by considering expansions on functions of infinite 
order of slow-varyness. As soon as we choose a finite order, as shown for ^ = 2, we may obtain 
nonlocal discrete equations as the integration of the equations defining the auxiliary harmonics 
may not be exact [28]. 

This work still leaves many open problems. Work is in progress on: 

• The proof of the Conjecture for some relevant classes of lattice equations. 

• The construction of higher orders in the multiscale expansion and on the definition of 
an order of integrability of a lattice equation (in the spirit of the paper [10]). 

It is worthwhile to notice that a meaningful classification of integrable Z^-lattice equations can 
be obtained only by going at least to the fourth order in the multiscale expansion (indeed, the 
third order is not enough to show the non- integrability of the discretized KdV equation ()16p ) 

In the present paper we have considered the modulation of weak plane waves solutions of the 
discrete systems. Work is also in progress on the modulation of constant solutions. In such a 
case the multiscale analysis should provide a KdV-type equation instead of a NLS one [26]. 
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